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Abstract 

A mean-field potential version of the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) 
model is used to investigate the production of strange baryons, especially the As and As, from heavy 
ion collisions at SPS energies. It is found that, with the consideration of both formed and pre- 
formed hadron potentials in UrQMD, the transverse mass and longitudinal rapidity distributions of 
experimental data of both As and As can be quantitatively explained fairly well. Our investigation 
also shows that both the production mechanism and the rescattering process of hadrons play 
important roles in the final yield of strange baryons. 
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The heavy ion collision (HIC) is the only way of the human being at present to explore 
the properties of nuclear matter at supranormal and subnormal densities and/or high tem- 
peratures. The properties of the sub-structure and the dynamics of the nucleus/nucleon 
could be even discovered by high-energy HICs which has being explored experimentally 
by the Schwerlonen Synchrotron (SIS) at Gesellschaft fur Schwerionenforschung (GSI, Ger- 
many), the Alternating Gradient Synchrotron (AGS) at the Brookhaven National Labora- 
tory (BNL,USA), the Super Proton Synchrotron (SPS) at the European Organization for 
Nuclear Research (CERN, Switzerland), the Relativistic Heavy Ion Collider (RHIC) at BNL, 
and the Large Hadron Collider (LHC) at CERN. Theoretically, the most important property 
of the nuclear matter relates closely to its equation of state (EoS) jl-4|. The stiffness of the 
5oS at both low and high densities have also being received much attention in past decades 
-I 111]. At high densities, the stiffness of the EoS will definitely decide the order and the 
level of the phase transition to the quark-gluon plasma (QGP) 12H16|. It is a consensus 
that any final conclusion of it should only be drawn by large numbers of comparison works 
between theoretical simulations and experimental measurements. 

In the past, several signals - such as charmonium suppression, relative strangeness en- 
hancement, energy loss of hard partons, etc - of the (phase) transition to the deconfined 



phase have been observed in HICs at SPS energies [17H22|. However, none of them gives 
100% undoubt conclusion of the probable phase transition. Especially, the strangeness pro- 
duction and enhancement is still a matter of argument. It is quite necessary to investigate 
the strange hadrons with a microscopic transport model which might give deeper insights 
into the whole process of the production and the transport of all hadrons. The Ultra- 



relativistic Quantum Molecular Dynamics (UrQMD) |l5l . |23M26| is one of the most suitable 
microscopic transport models which have been worked fairly successfully in this field for 
more than 10 years in addition to the Relativistic Quantum Molecular Dynamics (RQMD) 



the Parton-Hadron-String Dynamics (PHSD) |28j and A MultiPhase Transport model 
(AMPT) y models. 

However, based on previous comparison of UrQMD calculations with data, we know that 
there are some disagreement about the strange hadron production at AGS and SPS energies 
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311 ] . As for hyperons, first, the yields of hyperons are somewhat overestimated in 



UrQMD cascade calculations using versions earlier than 2.1. Starting from version 2.1, the 
UrQMD group considers additional high mass resonances which leads to a smaller yield of 
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hyperons so that a nice agreement with data from central Pb+Pb collisions at SPS energies 
was seen in previous calculations j^J. Second, the resulting mean transverse momenta of 
hyperons were found to be too low as compared with experimental data when using the 
version earlier than 2.1 [31]. It is interesting to see that this problem is also partly cured by 
the newest version 2.3. Third, even though the former two problems have been cured partly 
with the updating cascade version, the yield of anti-hyperons decreases at the same time 
which makes the comparison with data become even worse {31]. These problems attract 
our attention and will be discussed in this paper with a mean-field potential version of the 
UrQMD model. 

In this paper, after a brief introduction of the UrQMD and its recent updates, the pro- 
duction and the evolution of strange hadrons as well as anti-protons at SPS energies are in- 
vestigated with the mean-field potential version. For comparison, the corresponding cascade 
calculation results are also shown. The experimental data are taken from NA49 collaboration 
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341 ] . Finally, a conclusion and outlook is given. 



The UrQMD model is a microscopic transport approach based on the covariant propaga- 
tion of constituent (anti-)quarks and diquarks accompanied by mesons and baryons, as well 
as the corresponding anti-particles, i.e., full baryon-anti-baryon symmetry is included. It 
simulates multiple interactions of ingoing and newly produced particles, the excitation and 
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fragmentation of color strings and the formation and decay of hadronic resonances 
Besides the cascade mode in which all particles are treated to be free streaming between 
collisions, it is also necessary to incorporate the mean-field contribution for a complete dy- 



namic transport 15|, |35|, |36|. Since the UrQMD model inherits the basic treatment of the 
baryonic equation of motion in the QMD model {37], the consideration of the mean-field 
contribution is also logical. In order to properly describe the physical process in HICs at SIS 

-3- 
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With the 



energies this term has been treated carefully in the UrQMD model before 
increase of beam energy from SIS, AGS, SPS, up to RHIC, the dynamics of the transport 
has being attracted more and more attention due to the fact that quite a few of discrepancies 



between cascade calculations, of any microscopic transport models, and experimental 



have been shown in these beam energy regions, such as the collective flows 
nuclear stopping 
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data 



25 



40||, and the HBT two-particle correlation [36l. |42|. 



4l|, the 



Recent update of the UrQMD model follows three different routines. The first one is called 
as "cascade version" which is to modify the cascade process of UrQMD. This is the main 
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routine and the newest official version is v2.3 [3JJ. The second one is called as "mean-field 
potential version" which further considers the mean-field potentials of both formed and pre- 
formed hadrons [36] . The third one is called as "hybrid version" with which the microscopic 
transport process is incorporated with a macroscopic hydrodynamics 15|, |43|. So, by using 
the 'hybrid version' one can compare calculations with various EoS during the hydrodynamic 
evolution and those with the pure cascade calculations within the same framework. The first 
routine is of importance since it supplies a better basis as a cascade process. The latter two 
versions consider more deeply about the dynamic process of all particles but with different 
strategies: The second routine bases on the same structure of the equation of motion at low 
SIS beam energies, in which both the mean field and the collision terms should be taken into 
account explicitly. While the third routine replaces the UrQMD dynamic process in between 
the initialization and the hadron rescattering processes by a hydrodynamic one. Therefore, 
one sees clearly that, with the rapid development of researches on the dynamic process of 
high-energy HICs, the modifications on the transport model are still quite frequent and 
effective. 

n 

Starting from the version 2.0, similar to HSD, the PYTHIA |44| is considered in order 
to treat the initial hard collisions more carefully, which is important for HICs at high SPS 



and RHIC energies. In the version 2.1 25], similar to RQMD, the high-mass resonances 



are re-treated in order to give higher mean transverse momenta of most particles. The 



newest version 2.3 31] is then brought out with some other minor changes. The mean-field 



potential version is based on the cascade version v2.1 but not v2.3 which is partly due to 



the minor difference between them. As for t 
conventional potentials for formed hadrons 



re mean-field contribution, in addition to the 



35], the mean- field potentials for pre-formed 



36], i.e., the Yukawa, 



baryons from string fragmentation have also been taken into account 
the Coulomb, and the momentum dependent terms are neglected but the similar density 
dependent (Skyrme-like) term as the formed baryons is used, which reads as 

U(p h /p ) = a(^) + b(^y, (1) 
Po Po 

while a reduction factor (2/3) is considered for pre-formed mesons due to the quark- number 
difference. In Eq. (0Q) p Q = 0.16/mr 3 is the normal nuclear density, a, b, and g are param- 
eters, in this work for the SM-EoS, they are —110 MeV, 182 MeV, and 7/6, respectively. 
The ph is the hadronic density, which reads as 



Ph = y] QCjPij (2) 

where Cjj = 1 for baryons, 2/3 for pre-formed mesons, and for formed mesons, pij is a 
Gaussian in coordinate space. 

As in Ref . 45] , the relativistic effects on the relative distance r ^ = i - , — Xj and the relative 
momentum = p« — Pj employed in the two-body potentials (Lorentz transformation) are 
considered as follows: 

f 3 = 4+7^-/%) 2 , (3) 

2 9 
<YY) — Til 

rt=rt-(E t -E j f + 1 U 1 ^) 2 - (4) 
In Eqs. [3] and H] the velocity factor and the corresponding 7-factor of i and j particles are 



defined as = (p?. + Pj) /{Ei + Ej) and 7^ = 1/y 1 — Furthermore, a covariance-related 
reduction factor for potentials in the Hamiltonian, rrij/Ej, was introduced in the simplified 



version of RQMD model 



46| and adopted in this work as well. 



We have found from the HBT correlations of two-particles at AGS, SPS, and RHIC 



energies that the mean-field potentials of both formed and pre 



for explaining the sources of the HBT time- related puzzle 36l l47| . Furthermore, by using the 



brmed hadrons are essential 



hybrid model, it is found that the equation of state at high densities should be somewhat 



stiff in order to explain the HBT data at SPS energies [15j. Although, the stiffness of 
the EoS is still with some uncertainties, it is no doubt that the dynamic process deserves 
more investigations. We notice that the average hadronic density at a central zone (e.g., 
Rem. < 5/m) from central Pb+Pb reactions at SPS energies reaches 4-8 times normal 
nuclear density at the early stage of the collision. And, it is known that in the same beam 
energy region and at the early stage, the string process starts to replace the resonance-decay 
process for producing new particles in the UrQMD model description. Fig. [1] shows the 
time evolution of the ratio of pre-formed hadrons and all hadrons ("Rpre-formed/aii") in the 
zone R c . m . < 5fm from the central Pb+Pb reaction at 40A GeV (solid line) and 158A GeV 
(dashed line), respectively. It is seen clearly that the main production mechanism at t < 5 
fm/c is the string excitation and fragmentation. And it is up to t ~ 10 fm /c that this 
production mechanism still plays visible role. Hence, the pre-formed hadron potentials will 
definitely provide a large contribution to the early pressure. Accordingly, the rescattering 
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FIG. 1: Time evolution of the ratio of pre-formed hadrons and all hadrons ( ll R pre - formed/ all") i n 
the zone R c . m . < 5/m from central Pb+Pb reaction at 40^4 GeV (solid line) and 158A GeV (dashed 
line), respectively. 

of new produced particles will be influenced to some extent, which will be checked by the 
current investigation. 

Fig. [2] shows the transverse mass (m t — mo, where rrit = \Zpf~+ ano ^ m o is the mass 
of the particle at rest, p t = ^Jpi + Vy is the transverse momentum of the particle in the 
center-of-mass system) spectra of As (= A + S°, same as experimental data) at mid-rapidity 
(\y\ < 0.4) for central Pb+Pb reactions (a/a T < 5%) at AOA GeV (left plot) and at 158/1 



GeV (right plot). The data are taken from Ref. [34J for comparison. We find that the 
cascade calculation of the UrQMD version 2.1 already reproduces the data well within error 
bars except those at very low and high transverse masses. We also noticed that if the newest 
version 2.3 is used this situation becomes even worse at large transverse masses, which might 
be due to the minor changes of the double strange diquark suppression factor and the single 
strange diquark suppression factor. If we use the older version 1.3, as we have known before 



3jJ , the slope of the transverse mass spectra is steeper than the newest one. Hence the 



careful treatment of the high-mass resonances is very important. With the consideration of 



6 



c/ct t :<5% 




0.2 0.4 0.6 0.8 1.0 1.2 0.2 0.4 0.6 0.8 1.0 1.2 

m t -m (GeV) 



FIG. 2: Transverse mass spectra of As (= A + E ) at mid-rapidity (\y\ < 0.4) for central Pb+Pb 
reactions (<j/<jt < 5%) at AO A GeV (left plot) and at 158.A GeV (right plot). The experimental 



data are shown by stars (taken from Ref. 



341 ]) while calculations with mean- field potentials (dashed 



lines) are compared with those with a cascade mode (solid lines). 

both formed and pre-formed hadron mean-field potentials, it is found that the slope becomes 
more flat so that one can further describe the data at both the low and the high transverse 
masses. Therefore, the treatment of the mean field contributions is also of importance. 

Fig. [3] depicts centrality dependence of the rapidity distribution of A yields for Pb+Pb 
collisions at 40 A GeV (upper plots) and 158 A GeV (lower plots). Results within five central- 
ity bins from central (cr/crr < 5%) to semi-peripheral (33.5% < a/ar < 43.5%) collisions are 
shown from the right to the left. First of all, it is seen that the effect of potentials is stronger 
in more central collisions, which is clear. Second, the potential effect is mainly to suppress 
the yield of hy per ons at mid-rapidity which was also seen in previous calculations at lower 



beam energies 



301 ]. It implies again that the dynamic transport of the hyperon after its pro- 



duction is important. Third, the suppression effect on the yield of hyperons provides us with 
a better fit to data especially at the most central collisions. We also find that calculations 
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FIG. 3: Rapidity spectra of As for Pb+Pb collisions at 40^4 GeV (upper plots) and 158A GeV (lower 
plots) for 5 different centrality bins which are shown from the right to the left. The experimental 



341 ]. Calculations with mean-field potentials 



data in each plot is shown by stars and taken from 
and with a pure cascade mode are shown by lines. 

with the version 2.3 also present similar results to the ones with potentials. Finally, it is also 
seen that in the most central Pb+Pb collisions at 158A GeV, the rapidity distribution of A 
yields calculated with potentials is slightly expanded. This phenomenon should be related 
to the stronger early pressure introduced by the pre-formed hadron potentials and will be 
discussed later-on. 

Fig.HJshows the rapidity spectra of As (= A+£°) for central Pb+Pb collisions at AO A GeV 
(left) and 158 A GeV (right). Calculations with and without mean- field hadron potentials 



are compared to the data taken from [34]. Contrary to the rapidity distribution of As, it is 
interesting to see that the yield of As at mid-rapidity is driven up to fit the data quite well 
when the potentials are considered in calculations. 

In order to understand the nice fitting results of both A and A spectra in both longitu- 
dinal rapidity spectra and transverse mass spectra, and to check the fitting results of other 
particles, we show in Fig. [5] the rapidity spectra of S~s (left) and anti-protons (right). The 
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FIG. 4: Rapidity spectra of As for central Pb+Pb collisions at 40A GeV (left) and 158^4 GeV 
(right). Experimental data and the data reflected at mid-rapidity are shown by solid and open 
stars [34]. Calculations with and without the mean-field potential are shown by lines. 

central Pb+Pb collision at 40A GeV is chosen as an example since the experimental data 
are available for both particles. It is seen clearly from Fig. [5] that the calculated S~ yield 
is suppressed with the consideration of potentials, which is the same as As. While the cal- 
culated yield of anti-protons is enhanced which is the same as anti-As. At 158A GeV, the 
same potential effect is seen for both particles. 

Therefore, it is necessary to check the time evolution of these produced particles in order 
to understand the potential effects on both hyperons and anti-particles. Fig. [6] exhibits the 
calculated transverse mass spectra of H~s (left) and anti-protons (right) at two evolution 
time points t = 3 fm/c and 30 fm/c. The central Pb+Pb reaction at the beam energy 158A 
GeV is chosen in order to have more particles during the time evolution. From Fig. [6J we 
find that shortly after the time when the two Lorentz-contracted nuclei have passed through 
each other 48[ (e.g., ~ 3 fm/c and ~ 1.5 fm/c for the Pb+Pb system at 40A GeV and 158A 
GeV, respectively), here we set t = 3 fm/c, the potential effect on E~ production can be seen 
at high transverse mass, but it is more obvious in the anti-proton spectra. It should be due 
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FIG. 5: Rapidity spectra of E s (left) and anti 



GeV. Experimental data are shown by stars 
potential are shown by lines. 



protons (right) for central Pb+Pb collisions at 40^4 
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331 ] . Calculations with and without the mean- field 



to the fact that anti-particles are mainly produced at earlier collisions within a more dense 
region. Therefore, a higher early pressure introduced by the pre-formed hadron potential 
leads to stronger emission of anti-protons with high momenta. With the increasing time 
from t = 3 fm/c to t = 30 fm/c (when most of rescatterings of particles have ceased), we 
find that in the cascade mode there are still a large amount of new S~s produced, while the 
production of anti-protons in this period becomes much less important. The time evolution 
of total yields of particles can be seen more clearly from Fig. [7J in the cascade mode and with 
the time increasing from 3 fm/c to 30 fm/c, the total yield of S~ is doubled increased, while 
the yield of anti-protons is about 30% decreased due to the well-known strong annihilation 
effect. It is also found that after 12 fm/c the yield of anti-protons is almost saturated. With 
the consideration of formed and pre-formed hadron potentials, however, more anti-protons 
are survived which are mainly due to their higher momenta obtained at the early stage of 
reactions. If we switch off the pre-formed hadron potentials but keep the formed ones, it 
is seen from Fig. [7] that the time evolution of anti-protons is almost the same as that with 
the cascade mode. It implies that the mechanism of string excitation and fragmentation is 
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FIG. 6: Calculated transverse mass spectra of H _ s (left) and anti-protons (right) for central Pb+Pb 
collisions at 158A GeV and at two time points: 3 fm/c and 30 fm/c. Calculation results with and 
without the mean- field potential are shown by lines with symbols. 

essential to the production of anti-protons, which certainly is reasonable. For strange particle 
production, in addition to the string mechanism, the rescattering process of hadrons are also 
of importance, which can be understood clearly from the following features shown in Figs. E] 
and [7J (1) the rapid increase of the S~ yield during the time 3 — 30 fm/c, (2) the suppression 
effect of potentials on both the yield mainly at the low transverse masses shown in Fig. [6] 
and the total yield shown in Fig.[7]at t = 30 fm/c, and (3) the contribution of formed hadron 
potentials to yield (the line with half-open circle in Fig. [7j). 

Thus, up to now, we have understood more deeply the reason why the mean-field po- 
tentials contribute to the nice fitting results of both A and A spectra shown in Figs. EH 
However, it is also found from Figs. [5]{T] that the modification of the UrQMD model should 
carry on with more endeavors. First, the production of multi-strange baryons is still not 
satisfactory. The yield at mid-rapidity is underestimated at SPS energies in all versions of 
UrQMD calculations (including the newest version 2.3). This problem can not be solved 
solely by the change of the double strange diquark suppression factor and/or the single 
strange diquark suppression factor. Since the strange-hadron related cross sections are nor- 
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FIG. 7: Time evolution of total yields of H~s and anti-protons for central Pb+Pb collisions at 
158^1 GeV. Cascade calculation results (line with open symbols) are compared to the results with 
potentials of both formed and pre- formed particles (lines with solid symbols), and to the results 
with potentials of only formed particles (lines with half-open symbols). 

mally large, the medium effect of cross sections on the production and the rescattering of 
multi-strange baryons should show importance and deserves more investigations. Second, 



the production of (anti-) protons at AGS and SPS energies are still an open problem 3ll.l40|. 



Although the net-proton yield and the rapidity distributions of (anti-)protons are better de- 



scribed 



3 



40| by the mean-field potential version, a further refinement is still required. The 
dynamic process of particles in the dense medium induced by both the new (QGP) phase 
and the hadronic phase afterwards should be investigated with a theoretical breakthrough. 

To summarize, the production of strange baryons, especially the As and As, are investi- 
gated with the mean-field potential version of the UrQMD model for HICs at SPS energies. 
In this version of UrQMD, in addition to the formed hadrons, the mean-field potentials of 
pre-formed hadrons are also considered which are similar to those of formed hadrons. It is 
found that with the consideration of potentials, the transverse mass and longitudinal rapid- 
ity distributions of experimental data of both As and As can be quantitatively explained 
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fairly well. The early strongly repulsive and the later-on attractive forces introduced mainly 
by the pre-formed and the formed potentials respectively lead to more A emission but less A 
emission at freeze-out. Our investigation also shows that the hadronic rescattering process 
is still important for HICs at SPS energies. 
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